Single and low-lying states dominance in two-neutrino double-beta decay 
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A theoretical analysis of the single-state dominance hypothesis for the two-neutrino double-beta 
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decay rates is performed on the examples of the double-beta decays of Mo, Cd, and Te. 
We also test the validity of an extended low-lying-state dominance that takes into account the 
contributions of the low-lying excited states in the intermediate nucleus to the double-beta decay 
rates. This study has been accomplished for all the double-beta emitters for which we have experi- 
mental information on their half-lives. The theoretical framework is a proton-neutron quasiparticle 
random-phase approximation based on a deformed Skyrme Hartree-Fock mean field with pairing cor- 
^ , relations. Our calculations indicate that there are not clear evidences for single- or low-lying-state 

O ' dominance in the two-neutrino double-beta decay. Finally, we investigate the single electron energy 

' distributions of the outgoing electrons in the double-beta decay processes with an exact treatment 

of the energy denominators, which could help to a more comprehensive analysis of NEMO-3 data. 

m ■ 

PACS numbers: 21.60. Jz;23.40.Hc;27.60.+j 

I. INTRODUCTION 

o 

(-H I The nonzero mass of the neutrino has been recently confirmed by neutrino oscillation experiments where 
neutrinos created in a flavor eigenstate are subsequently found to be in various flavors. However, these experiments 
can only probe differences of squared neutrino masses, but not the absolute neutrino mass scale, which remains stiU 
unknown. The neutrinoless double-beta (0i//3/?) decay nuclear process 0, d, 01 is considered to be one of the most 
Qv^ suitable candidates to provide this information. decay is a lepton-number violating mode, which is forbidden 

^-H in the Standard Model and that only occurs if neutrinos are massive Majorana particles. The half-life of this process 
CO involves an effective neutrino mass that could be eventually extracted from the measured half-lives. 

' Since the double-beta decay process is necessarily dependent on the nuclear structure properties, lack of accuracy in 
. the determination of the nuclear matrix elements involved in the process is a source of uncertainty in the information 
T-H on neutrino properties that can be extracted from Qv(3l3 decay experiments. Contrary to the Qv[3l3 decay, the doublc- 
00 . beta decay with the emission of two (anti)neutrinos {2vf3(3) can proceed as a perturbative process in the Standard 
■ Model and has indeed been observed in several nuclides. This Standard Model allowed form of the decay is used as 
J> ] a test for the nuclear models considered. The operators involved in the Ov and 2v modes are different. In the 2v 
decay mode the operator connects the initial and final 0'*' states via virtual 1^ transitions to the intermediate nucleus. 
In the Qv decay mode the transitions to the intermediate nucleus may take place to many different multipolarities. 
^ i Nevertheless, the underlying nuclear structure involved in both processes is similar. Success in describing the 2v decay 
" " ' mode is a requirement for a reliable calculation of the nuclear matrix elements related to the Qvfip decay. 

The nuclear structure calculation involved in the 2i>(3P decay is not an easy task. Being a second order process in 
the weak interaction, the nuclear matrix elements involve a summation over the full set of virtual nuclear 1+ states 
in the odd-odd intermediate nucleus. Two leading microscopic nuclear models are commonly used to evaluate the 
nuclear matrix elements involved in the decay process. They are the nuclear shell model [5, 6] and the proton-neutron 
quasiparticle random-phase approximation (pnQRPA) 0, H, • The shell model approach takes into account all 

possible correlations but within a restricted valence space. It successfully describes the low-lying excited states, but 
has difficulties with the description of states at high excitation energies, in particular, in the region of the Gamow- 
Teller (GT) resonance for open-shell medium and heavy nuclei, where most of the double-beta emitters are located. 
On the other hand, pnQRPA calculations do not have these problems related to the limited model space, but relevant 
ground state correlations might be missing. In addition, these calculations are very sensitive to model parameters. 
One problem of the theoretical 2^(3/3 decay studies is whether the contributions of higher-lying states to the 2i//3/3 
decay amplitude, which are apparently disfavored by large energy denominators, play an important role. As we have 
mentioned, this has important consequences on the reliability of the shell model and pnQRPA to describe the nuclear 
matrix elements. 
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In the eighties, Abad and collaborators suggested ^ll'| that in those cases where the ground state of the odd-odd 
intermediate nucleus is a 1"*" state reachable by a GT transition, the transition through it could account for the major 
part of the 21/(3(3 matrix element. This conjecture is called the single-state dominance (SSD) hypothesis. If SSD 
hypothesis is confirmed, the half-lives for 2zy/3/3 decay could be determined from single (3~ and electron capture (EC) 
measurements, experiments which are in principle much easier to carry out than the elaborated and time-consuming 
double-beta decay experiments. On the theoretical side, the possible realization of the SSD hypothesis for the ground- 
state to ground-state transitions would lead to a drastic simplification in the theoretical description of the intermediate 
nucleus, since now only the lowest 1"*" wave function has to be calculated. If SSD is valid, this would increase the 
reliability of the small valence spaces used in the full shell model calculations. 

From the experimental point of view, early tests of the SSD hypothesis were evaluated on the 2i//3/3 decay of ^^^Mo 
and ^^^Cd, two cases where the ground state of the intermediate odd-odd nucleus is J'^ — 1+. In Refs. and 
[l^ . the EC half-lives of the ground state of the intermediate nucleus, ^'^"Tc and -'^^^In respectively, were measured. 
Combining these measurements with the already known (3~ half-lives, the SSD hypothesis was confirmed to some 
extent on those examples. Experimental information on the EC and (3~ decay of ^^^I is also available 

Theoretically, the SSD hypothesis was systematically studied in Ref. [l^ within a pnQRPA formalism. The analysis 
performed there showed that the SSD hypothesis is realized in most cases but the mechanism leading to this property 
is not unique. In some instances the SSD is realized through a true dominance of the first 1+ virtual state, while in 
other cases it is realized due to cancellations among the contributions from higher l ying 1"*" states of the intermediate 
nucleus. However, one should notice that the theoretical formalism used in Ref. [l5|, is complemented with some 
adjusting of the single-particle levels around the proton and neutron Fermi surfaces to reproduce the observed single- 
quasiparticle spectrum. In addition, the scaling factor of the strength in the particle-particle residual force is fixed 
in each nucleus by optimizing the agreement with the experimentally known EC and (3~ decay rates. This indicates 
that the theoretical evidences of the SSD hypothesis claimed in Ref. [T^ rely on the experimental ones. 

In Ref. [iBl the SSD hypothesis was studied considering an exact treatment of the energy denominators of the 
perturbation theory. It was found that the 2j//3/3 half-lives are reduced by a factor of about 20% for ground-state 
to ground-state transitions. It was also shown that by measuring the single-electron spectra and/or the angular 
distributions of the emitted electrons, the SSD hypothesis could be confirmed or ruled out experimentally in the near 
future. 

In the last years, experimental information on the low-lying Gamow- Teller strength coming from charge-exchange 
reactions has been collected for several nuclei relevant to double-beta decay processes [17, 18, 19, 20, 21] ■ In particular, 
(p, n) reactions were used in Ref. [iH to study the transitions ^^*Te(0''") — > ^^^I(l+), the {^¥Le,t) reaction was used in 
Ref. [lil to study the transitions i°°Mo(0+) -> io''Tc(l+) and "6Cd(0+) ^^Hn{l+). The (p, n) reaction was also 
used in Ref. ^ to study "SCd(0+) ^^^ln{l+). Similarly, the reaction {d,'^Re) was used in Ref. [IH to study the 
transitions ^^^Sn(0"'') ^^®In(l+). These measurements complement the experimental information already available 
on the ft values for the transitions from the l"*" ground state of the odd-odd intermediate nucleus to the 0+ ground 
states of parent and daughter nuclei. In addition, they give information on the transitions to the l"*" excited states of 
the intermediate nucleus. 

In this paper we use a formalism based on a pnQRPA approach with a selfconsistent quasiparticle basis obtained 
from deformed Skyrme Hartre-Fock calculations with pairing correlations [2^ . On top of that we include separable 
residual spin-isospin interactions in both particle-hole and particle-particle channels [23|. Within this formalism we 
calculate the GT strength distributions of the single beta branches as well as the 2v(3P matrix elements corresponding 
to the double-beta emitters ^°°Mo, ^^^Cd, and ^^^Te. With these results, we first evaluate to what extent the SSD 
hypothesis or the more relaxed low-lying states dominance (LLSD) hypothesis are fulfilled within our theoretical 
approach. Comparison with the available experimental information on both single-beta branches, EC / (3^ , and 2v(3(3 
decay is also made. Finally, we extend the LLSD analysis to the rest of the observed 2v(3(3 emitters, where the ground 
state of the intermediate nucleus is not 1^, and we evaluate the contributions to the 2v(3(3 half-lives considering an 
increasing range of excitation energy in the intermediate nucleus. Former calculations of running sums of matrix 
elements as a function of the excitation energy of the 1+ intermediate nucleus were presented in Refs. [2^. [25]. In 
Ref. [13] the results presented were for ^^Ge and were obtained within a phenomenological Woods-Saxon mean-field 
approach. In Ref. [25| results were presented for ^^Ar, '^'^Fe, and ^*Ni within a shell model approach. 

The paper is organized as follows. In Sec. II we present briefly the main formalism and discuss various approxima- 
tions. In Sec. Ill we compare our results with experiment in the cases of ^"°Mo, ^^^Cd, and ^^^Te and analyze the 
possible realizations of SSD and LLSD hypotheses. We also study in this section the single electron energy distribu- 
tions of the outgoing electrons in the 2v(3(3 decay processes with an exact treatment of the energy denominators. In 
Sec. IV we study the contributions of the different energy ranges to the 2v(3(3 matrix elements in all the observed 2v(3(3 
emitters and discuss them in terms of LLSD. Finally we present in Sec. V the summary and the main conclusions. 
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II. BRIEF DESCRIPTION OF THEORETICAL CALCULATIONS 

The double-beta decay is a nuclear process characterized by a change of two units in the nuclear charge, while the 
atomic mass number remains unchanged. Here we focus on the 2v(3~ I3~ decay mode where two antineutrinos and 
two electrons are emitted, 



(Z, A)^+ ^ (Z + 2, + 2e- + 2Pe 



(1) 



Only the basic expressions are written here, more details can be found in Refs. [2J,|26|. The half-life of the 21/(513 
decay can be written as 
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where G"^"^^ is the phase-space integral [5, ^ . The nuclear matrix element M^r^^ contains all the information on the 
nuclear structure involved in the process. The theoretical formalism used here is the pnQRPA based on a selfconsistent 
mean field that allows for pairing and deformation. We refer to this formalism in short as deformed HF-|-BCS-|-QRPA. 
In this formalism the nuclear matrix element can be written as, 
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where the intermediate states are the 1"*" excitations of the initial and final nuclei in the laboratory frame. We use 
the Bohr-Mottelson factorization [131 for the total angular momentum states /+ in the laboratory frame. The index 
fJ-iip-f) Eq. ^ contains the Ki{Kf) third component of the angular momentum of the QRPA state in the intrinsic 
system. 
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(4) 



with the intrinsic state 



Ki 



(5) 



with 100, i) the QRPA vacuum of the initial state and F^' the phonon QRPA operators. \(j>^') are the time reversed 
of |(/>^')- Similarly, the 0+ initial and final states are of the Bohr-Mottelson form 



|0+0) = -^Dl°{mo,^) , 



(6) 



and similarly for the final state /. Integrating over the angular fl variables, we finally obtain 



M, 



GT 
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(7) 



As it is well known (see Ref. [28'] and references therein) the Bohr-Mottelson factorization of the wave function gives 
the adiabatic limit of the exact angular momentum projection and is a very good approximation for well deformed 
nuclei. First order corrections to the Bohr-Mottelson factorization due to angular momentum projection are typically 
less than !/{,]']), where (J^) is larger than 100 for well deformed nuclei [2^. In Eq. ([7]) u;^'{u;^') are the RPA 
excitation energies of the intermediate l"*" virtual states with respect to the initial (final) nucleus, with energy Ei{Ef ). 
The energy denominator is the average excitation energy of a couple of these intermediate states. The indices mi, 
ruf run over all the 1"^ states of the intermediate nucleus {Z + 1, A) in the transition, and K = 0, ±1 are the 
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possible angular momentum projections of these intermediate states. The overlaps are needed to take into account 
the nonorthogonality of the intermediate states reached from initial and final ground states. Their expressions can be 
found in Refs. 

To compute the matrix element wc proceed as follows. Each initial and final even-even ground state is described as a 
pnQRPA vacuum. For each [N, Z) nucleus, a BCS ground state is independently obtained from a selfconsistent mean- 
field calculation on a large axially-symmetric harmonic oscillator basis, using a density-dependent effective interaction 
and a phenomenological pairing-gap parameter. The nucleon-nucleon effective force SLy4 [sol] is used to generate the 
HF mean field, and the BCS equations are solved at each iteration with a value of the pairing gap obtained from 
even-odd experimental mass differences between neighboring nuclei [3l| . When the iterative process converges, a set 
of single-particle states is obtained with their corresponding wave functions, energies, and occupation probabilities. 
It is worth noticing that the quadrupole deformation of the nucleus in its ground state is obtained selfconsistently 
as the shape that minimizes the nuclear energy. Separable Gamow- Teller particle-hole (ph) and particle-particle (pp) 
residual interactions are included and treated in pnQRPA [23| . The corresponding coupling constants Xph and Kpp 
have been used in this work with values Xph = 20 /A MeV and Kpp — 3/A MeV, which give a good description of GT 
properties when SLy4 Skyrme force is used, as we shall show in next section. 

Due to the many-body approximations performed in the calculation, the intermediate l"*" states computed from 
the initial and from the final ground states are different, and a projection from one basis to the other is in order. 
To get the 2i^/3P decay matrix element, every possible transition through all the intermediate 1+ states is taken into 
account. The amplitude of each transition is divided by an energy denominator, as shown in Eq. ([7]). The energy 
denominators written in this way are obtained by replacing the lepton energies E^i + E^j, where E^ {E^) is the energy 
of the emitted electron (antineutrino) , by an average quantity 0, H, 0] 



E^i+E^jK-Qpp + m^ = -{E,-Ef), (8) 

where Qpp is the Q-value of the double-beta decay [Qpp = Ei — Ef — 2me). Ei and Ef are the ground-state energies 
of parent and daughter nuclei, respectively. The consequences of this approximation on the SSD tests have been 
studied in Ref. We shall come back to this point in the next section. Different ways of writing these energy 

denominators, which are all equivalent, can be found in the literature. In particular one finds 



+ = + E,, - E, + = E^ - ^{E, + Ef) , (9) 

where Em is the energy of the intermediate nucleus in the state m. In the cases where the ground state of the 
intermediate nucleus is a 1"*" state, the denominator corresponding to the decay passing through this state can also be 
written as {Qp- + Qec)/2, where Qp- and Qec are the Q— values corresponding to its /3~ and EC decay, respectively. 

In this paper we study the contributions to the 2i//3/3 nuclear matrix elements from consecutive terms in Eq. ([7]) 
with increasing energy denominators. The main purpose is to check the validity of the SSD hypothesis by considering 
only the lowest-energy term in the summation in the cases of ^'^'^Mo, ^^®Cd, and ^^'^Te, which are emitters with 1+ 
intermediate ground states. However, it is also interesting to check the extent to which the contributions coming 
from low-lying virtual excitations are able to reproduce the total 2vj3(3 nuclear matrix elements in a more general 
way. Actually, the strong fragmentation of the GT strength within the deformed formalism suggests that a meaningful 
comparison of theoretical results with experiment should be done for the accumulated strength in a given energy range 
rather than a direct comparison of individual excitations. Thus, we analyze the summation in Eq. ([7]) considering 
only the contribution of the intermediate ground state (to test SSD), then the contribution corresponding to all the 
low- lying states (to test LLSD), and finally taking into account all possible contributions (to get the complete result). 

As mentioned before, in the cases where the ground state of the intermediate nucleus is a 1+ state, we have 
experimental information on the j3~ and EC decays of the ground state intermediate nucleus into the 0"*^ ground 
states of the daughter and parent nuclei, respectively. We also have complementary experimental information based 
on charge-exchange reactions that can be evaluated and compared to the 2i> f3P experimental nuclear matrix element, 
as well as to the theoretical calculations. The contribution involving low-lying states can also be evaluated from the 
experimental information including now charge- exchange reactions [l9l . [2l| , and can be compared again with the total 
2v(3f3 experimental nuclear matrix element, as well as to the theoretical calculations. 
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III. RESULTS FOR 2z//?/3 DECAYS WITH 1+ GROUND STATE INTERMEDIATE NUCLEI 



A. Single-state dominance 



We start this section by summarizing the experimental situation concerning the vaHdity of the SSD hypothesis in 
the 2v/3(3 decay of ^"^Mo, ^^^Cd, and ^^'^Te. The analysis is based on information extracted from ground-state to 
ground-state decay measurements of the intermediate nuclei ^°°Tc, -'^^^In, and as well as on information extracted 
from charge- exchange reactions. Table I, II, and III contain the GT matrix elements for ground state to ground state 
transitions in the cases A = 100, A — 116, and A = 128, respectively. The notation used corresponds to the reaction 
point of view, where the B{GT~) strength corresponds to the {Z, A) (Z + l, A) transition and the B{GT~^) strength 
corresponds to the {Z + 2, A) ^ {Z + 1, A) transition. Thus, in the case of the decay of the intermediate nucleus, we 
still call B(GT^) to the strength corresponding to the EC / (3^ decay, and B{GT'^) to the strength corresponding to 
the (3~ decay, according to the following scheme: 
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In this paper, results for B{GT) are given in units in which the neutron decay has B{GT)=3. The matrix elements 
are extracted from the cross sections in the case of charge-exchange reactions [2l| and from the log ft values in the 
case of the f3 decay : 



M{GT) 



3D ' 



1/2 



where D = 6147 and gA — 1.25 is the axial- vector coupling strength, 
within the SSD hypothesis is obtained from 



(10) 



The GT matrix element for the 2^(3/3 decay 



M'-'^^(SSD) M{GT-)M{GT+) ^ _l I 6D 



where the overlap in Eq. Q has been approximated by one. 

Table I shows the experimental Gamow- Teller matrix elements M{GT~) and M(GT+) for the transitions connecting 
the l"*" ground state in ^^^Tc with the 0+ ground states in ^°°Mo and ^°°Ru. The M(GT~) matrix elements were 



extracted from (•^He,t) charge-exchange reactions [19[ and from the ft^c value as given in Ref. [I4I. The M(GT^) 



matrix element was obtained from the ftp- value [30| in the decay ^"°Tc -'^^"Ru. The matrix elements M^j?'^ are 
obtained from the two possible combinations of data. The half-lives T^^^ are calculated using approximated (SSDl) 
and exact (SSD2) energy denominators (we keep the notation used in previous references [l6|). The experimental half- 
life is T^J^^ = 7.1 X 10^8 y [s^, which gives rise to a matrix element {Mq^^)^^^ = 0.241 MeV-\ using gA = 1.25. 
One should notice that the half-lives are independent on the gA value. As one can see in Table I, the agreement 
with the experimental half-life improves with SSD2 approach and it is particularly good when the matrix elements 
extracted from log ft are used. 

Table II is similar to Table I but for A=116. Experimental values for M{GT~) are taken from (^IIe,t) 19] and 
(p, n) charge- exchange reactions, and from the /isc value for the decay of ^^^In into the ground state of ^^^Cd 
reported in Ref. [l^. Experimental values for M(GT+) are from Ref. [2l|, where the GT strength was extracted from 
the ((i,^He) charge-exchange reaction ^^^Sn ^^^In. In this reference, the strength was normalized to recover at zero 
excitation energy the /3~ decay data ■ The results for Mq^^^ are again compared to the values extracted from the 
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measured 2z^/3/? decay half-life [111 {T^^^ = 3.0 x 10-^^ y, and a corresponding matrix element {M^^^)cxp = 0.127 

MeV~^, using gA = 1.25). Except for the first row that contains the values extracted from {^lie,t) and that clearly 
overestimates the experimental half-life by almost one order of magnitude, we observe that both SSDl and SSD2 
half-lives underestimate the measured value. 

Table III contains similar information for the case A = 128. Experimental values for M{GT~) are taken from the 
/^EC value for the decay of into the ground state of ^^^Te reported in Ref. [l3|, as well as from {p, n) reactions fl8l |. 
Experimental values for M{GT^) are taken from the ftp- value for the decay of ^^^I into the ground state of ^^^Xe 
[ll|. In this case both SSDl and SSD2 predictions overestimate the experimental half- life [sS] {T^^'^ = 2.5 x lO^** y, 
and a corresponding matrix element (Mq^^)cxp — 0.043 MeV~^, using qa — 1.25). 

B. Low- lying states dominance 

In Tables IV, V, and VI, for the cases A — 100, A = 116, and A — 128 respectively, we show theoretical results 
and experimental data corresponding to the B{GT^) strengths to the low- lying virtual states of the corresponding 
intermediate nuclei, as well as the contributions to 2i/f3/3 decay matrix element of the virtual transitions through them. 
A quenching factor gA,cS — O.SgA.hcc has been introduced in all the theoretical calculations reported in this paper, 
except for those shown later on in Fig. [SI In Table IV we do not have information on the low- lying B{GT^) strengths 
from charge-exchange reactions, only the log ft value is available. Thus, we only compute theoretical values for the 
2vl3(3 decay matrix element. We can see that the theoretical B{GT~) strengths accumulated up to 1.4 and 2.6 MeV, 
agree with the measurements [l^ within a 80-90%. On the other hand, the 2vl3l3 decay matrix element calculated up 
to 2.6 MeV of excitation energy in the intermediate nucleus accounts only for 65% of the experimental value. 

In the case of A = 116 we can combine the experimental information available from (•^He,f) ^] and from (c?,^IIe) 
[2ll | to evaluate the 2vl3l3 decay matrix element M'^r^^ up to 3 MeV, and we can compare the obtained values with 
the total matrix element evaluated from the measured 2v(3j3 half-life. Thus, we can test to what extent LLSD is valid 
in this case both theoretically and experimentally. First, we observe in Table V that the two single branches are 
reasonably reproduced by the calculations. In the case of i?(GT+) the agreement is apparent. In the case of B{GT~) 
the comparison to experiment is not so meaningful because there is no agreement between experimental data from 
different groups p^ . |20| . In particular, a very different value has been reported for the transition to the ground 
state. If we consider the new revised values given in Ref. [2^ (see Table II), the agreement with the calculations 
clearly improves. The 2vf3l3 decay matrix elements evaluated at increasing energies show first that they approach to 
the experimental value obtained from the 2vl3(3 half-life (0.127 MeV^^). We also observe that the theoretical results 
increase with energy to yield about a half of the experimental value at 3 MeV. Agreement improves when using the 
value of Ref. [l^]. Table VI contains the results for A = 128. As in the case of ^4 = 100, we do not have information 
on the low-lying B{GT^) strengths. Thus, we only compute theoretical values for the 2v(3l3 decay matrix element. 

Fig.[T]shows the accumulated single B{GT~) strength from -'^°°Mo — > -'^'"'Tc (upper panel), the accumulated single 
B{GT+) strength from lO'^Ru ^^^Tc (middle panel) and the 2vl3(3 GT matrix element from '^^^Mo ^""Ru (lower 
panel) in the low excitation energy range (up to 5 MeV). Available experimental data are shown by thick lines. 

As one can see in Fig. [T]the theoretical B{GT~) strength below 3 MeV is rather fragmented, but the total accu- 
mulated strength agrees well with experiment. Although not shown in the figure, it is worth noticing that the GT 
resonance measured in Ref. [l^ at 13-14 MeV with a strength B{GT^) ~ 23 is well reproduced in our calculation, 
yielding a peak centered at 16 MeV with a similar strength. The strength obtained from the log ft value also agrees 
with the strength calculated up to about 0.5 MeV. In the case of the B{GT^) strength the only information comes 
from the log ft value at zero excitation energy. The corresponding strength is again well reproduced by the calculation 
of the accumulated strength below 0.5 MeV. 

In the lower panel we show the theoretical calculation for the 2v[3l3 matrix element as a function of the excitation 
energy taken into account in the sum in Eq. ([7]). It is compared with the experimental matrix element extracted 
from the measured 2v[3(3 half-life. We compare our results with the values extracted from the measured log ft and 
charge-exchange reactions to check the validity of the SSD hypothesis graphically. We also show the total theoretical 
value when all the intermediate 1"*" states are considered in the calculation (see also Fig. |6]). As can be seen in the 
lower panel, a rapid increase in the 2v[3(3 decay matrix element takes place up to 2 MeV of intermediate excitation 
energy. It accounts for about 60% of the total matrix element. Eventually (see Fig. ^ the calculations reach the 
experiment when one considers the whole energy range. From the single-beta strengths one can see that there are no 
especially strong transitions to any low-energy intermediate state (in particular, to the ground state), which prevents 
the 2vf3j3 matrix element from showing the same effect. 

bummarizmg, m the 2vj3l3 decay of i°°Mo, the SSD hypothesis is fulfilled experimentally within 90%. The total 
theoretical 2z//3/3 matrix element agrees well with experiment, while LLSD calculations up to 2 MeV account for 60% 
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of the experimental 2i'(3j3 matrix element. The remaining 40% comes from contributions of states at higher excitation 
energy. 

Same results are shown for the transitions ^^^Cd -'^^^In (upper panel), ^^^Sn — > ^^^In (middle panel) and ^^^Cd 
^^^Sn (lower panel) in Fig. [21 In the upper panel we show experimental data for B{GT^) from (^He,i), together 
with the data at zero excitation energy extracted from (p, n), and from the log ft values. We also show the theoretical 
calculations that produce an increasing strength up to 1 MeV, which is very fragmented, but that finally reaches the 
most reliable experimental strength obtained from the measured log ft. In this case the GT resonance measured in 
Ref. [19] amounts to B{GT^) — 26 and it is located at 14.5 MeV. Our calculation yields about the same strength 
centered at 16 MeV. The experimental data for B{GT^) are from ((i,^He) charge-exchange reactions I2l|, where the 
value extracted from f3~ decay was used to normalize the strength. 

The B{GT^) strength in the middle panel shows that the calculation reproduces the total experimentally measured 
strength below 3 MeV. The distribution of this strength is more fragmented in the calculations but the general trend 
is similar. In the lower panel we compare again our theoretical results for the 2v[3l3 decay matrix element with 
experiment. In this case we compare with the experimental matrix element extracted from the measured half-life, as 
well as with the values extracted from the measured log ft values, and with the values extracted from charge- exchange 
reactions. We also show the total theoretical value (see also Fig. [S]), which is very close to the experiment. 

Theoretically, as in the case oi A — 100, a strict SSD or LLSD is not observed in the lower panel. The few 
contributions below 1 MeV of excitation energy account for a 50% of the final value of the double-beta decay matrix 
element. Experimental accumulated single strengths are roughly in agreement with the results shown, but the double- 
beta experimental data for the low-energy contributions are, as in the case oi A = 100, larger than our theoretical 
results. 

Summarizing, in the 2v(3j3 decay of ^^^Cd there is no clear evidence that the SSD hypothesis is fulfilled experimen- 
tally because of the spread of data coming from different sources. The experimental LLSD up to 2.5 MeV accounts 
for the experimental 2u(3[3 decay matrix element, while the calculations up to this energy account only for half of it. 
Since the total calculated matrix elements agrees well with experiment, the remaining contribution comes again from 
excited states at higher energies. 

Fig. [3] is the corresponding figure for A — 128. In the upper panel we show experimental data for B{GT^) from 
{p,n), together with the data at zero excitation energy extracted from the log/i values. We also show the theoretical 
pnQRPA calculations. In this case the GT resonance measured in Ref. [18] contains a strength B{GT~) = 34.24 
located at 13.14 MeV. Our calculation yields B(GT~) = 35 centered at 15 MeV. In the middle panel the experimental 
value at zero energy is obtained from log ft and agrees well with the calculation (note the different scale as 
compared to the upper panel). The lower panel contains the results for the 2vl3(3 decay matrix element. The meaning 
is the same as in the previous figures. We see that the total theoretical value (see also Fig. [5]), is very close to the 
experiment and it is not reproduced by SSD or LLSD hypotheses. 

A graphical representation of the calculations involved in Eq. ([7]) is shown in Fig. 01 featuring the theoretical l"*" 
spectrum of the intermediate nucleus ^°°Tc in the decay ^""Mo^^^^Ru in the left panel, the intermediate nucleus 
^^^In in the decay ^^^Cd^^^^Sn in the middle panel, and the intermediate nucleus ^^^I in the decay ^^^Te^^^^Xe in 
the right panel. The vertical scale is the excitation energy of the intermediate nucleus, the horizontal scale represents 
the B{GT) strength of the 1+ QRPA states of the intermediate nucleus as obtained from the GT~ transitions in the 
initial double-beta emitter (shown from the central axis to the left), and from the GT"*" transitions in the final double- 
beta partner (shown from the central axis to the right). The larger this strength, the more relevant the corresponding 
contribution to the total double-beta decay matrix element, provided there are states coming from the other branch 
of the decay with similar wave functions (giving rise to non-negligible overlaps) and carrying non- negligible strengths. 



C. Single electron energy distributions 

In Ref. [l^, the SSD hypothesis was tested using an exact treatment of the energy denominator in the nuclear 
matrix element, instead of the usual treatment of approximating the lepton energy by an average energy. It was 
shown that using the exact treatment of the energy denominator, the 2u(3j3 decay half-life in the case of ^'^'^Mo, is 
reduced by a factor of 20% for ground state to ground state decay. The single electron energy distributions of the 
outgoingelectrons were also analyzed, and the predictions of different assumptions, SSD and higher state dominance 
(HSD)|1Q], were compared. In this work we extend this analysis to the case of pnQRPA theoretical calculations of 
the 2v(3l3 decay rates with a correct treatment of the energy denominators and compare these predictions with those 
of SSD and HSD. 

The formalism needed to calculate the differential 2^(3(3 decay rates was developed in Ref. [13] . Here we only write 
the basic expressions used for the evaluation of the normalized decay rates. 
The differential 2uf3P decay rate to 0"*" ground state can be written as 
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dW = a2.F{Zf, E,i)F{Zf,E,2) M dil, 



(12) 



where — (G'/3.gyi)^mg/(647r'') and G/3 — GpcosOc [Gp is Fermi constant, 9c is Cabbibo angle). F{Zf,Ee) denotes 
the relativistic Coulomb factor. The phase space factor is given by 



TjEelPel Ee2Pe2 E^l £'^2 ^{Eel + Ee2 + Eyi + Ey2 + Ef — Ei) X 



dEel dEe2 dE^i dE^2- 
M. consists of the products of nuclear matrix elements: 



M = 



\JC + C\ 



;\}C-C\ 



where K.{C) denotes nuclear matrix elements with an energy denominator given by Km {Lr, 



(13) 



(14) 



Km = [Em — Ei + Eel + -B^l] ^ + [Em — Ei + Ee2 + ^ , 

Lm = [Em — Ei + Ee2 + E^i] ^ + [Em — Ei + Eel + E^2\ ^ ■ 



The normalized total decay rate can be written as 

P 

where the full decay probability W is given by 



W 
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1 dW 
WdK ' 



Ei-Ef-r. 



Ei-Ef-E^i 



F{Zf,Pe2)Pe2Ee2 dEe2 



F(Zf,Eel)pelEeldEel X 

Ei — E f — Ei=i — -Eg 2 

M El^El,dE,i . 



(15) 



(16) 



(17) 



In Fig. Owe present the single electron spectrum of the emitted electrons calculated within SSD, HSD, and total 
pnQRPA for the 2i//3/3 of ^°"Mo (upper panel) and ^^^Cd (lower panel). The different behavior observed at small 
electron energies for different assumptions is of the order of few percent. This accuracy is now accessible in the 
NEMO-3 experiment [s^, where a large amount of events for the (3(3 decay of ^"^Mo has been collected. 



IV. TEST OF LLSD IN DOUBLE-BETA EMITTERS 



As an extension of the analysis of the SSD hypothesis, we have also studied the contributions to the 2i/(3(3 decay 
matrix element through the low-lying intermediate states for all confirmed 2i//3/3 partners, also in the cases in which 
the ground state of the corresponding intermediate nucleus is not a 1+ state. To this end, we plot in Fig. [5] for all 
the nuclei under study the 2v(3(3 decay matrix elements as a function of the excitation energy of the intermediate 
nucleus taken into account in the calculation. As in the previous section, a deformed Skyrme-SLy4 HF mean field 
has been used to describe the initial and final ground states, with pairing correlations treated in BCS approximation. 
In general we use the ground states for initial and final nuclei with a shape that minimize the HF+BCS energy. The 
corresponding results are shown by the solid lines. In some instances like ^^Ti and ^^Se, the minima correspond to 
a spherical shape (/3 = 0), while experimental data |3Q] indicate that these nuclei have a non negligible deformation 
in the ground state {(3 ^ 0.15). In these cases we also show results obtained with the pnQRPA calculations based on 
the HF+BCS solutions corresponding to these non-zero deformations {(3 = 0.17 for '^^Ti and (3 — 0.14 for ^^Se). The 
corresponding results shown by dashed lines in Fig. [H] are seen to agree much better with experiment. It is important 
to remark that particularly in the case of ^*Ti, the HF-I-BCS energy has a rather shallow minimum in the deformation 
range from (3 = Q io (3 = 0.17, so that both solutions are equally plausible ground states of "^^Ti. The grey area in 
the plots indicates the experimental range of the 2v(3(3 decay matrix element. The two limiting horizontal lines have 
been deduced from experimental half-lives [ssj using two different values of the constant qa, namely 5^=1-25 (bare 
value) and (7^1=1. 00 (quenched value). In this case the theoretical results are not quenched. 
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Typically, the 2://?/? decay matrix elements shown in the plots increase with the excitation energy of the intermediate 
states until they reach a constant value right after a rather prominent final step located between 15 and 20 MeV. 
This last step, especially important for the double-beta partners with A = 128, 130 and 136, is related to the GT~ 
giant resonance. The final value of the 2v(3(3 decay matrix element lies within the experimental region in most cases. 
Important deviations are only found for the nuclei with A — 130 and 136. In the cases of ^ = 48 and A — 76, as 
already mentioned, the prolate shapes of the daughter nuclei give a better result than the spherical shapes. Concerning 
the low energy contributions to the matrix element, the fastest increase appears in the double-beta partners with A 
= 48, 96, 100, and 116, where around a 60% of the total matrix element is reached within an excitation energy range 
of 2 MeV, and to a lesser extent the A=150 case. In the other double-beta partners the contributions to the matrix 
element are more spread and the increase is slower. In general, important contributions appear from relatively high 
energies around the position of the GT resonance. Thus, although neither SSD nor LLSD hypotheses are clearly 
fulfilled in our theoretical results, one can see a tendency in most cases to exhaust at least 50% of the total 2vl3(3 
matrix elements in a low excitation energy region of 5 MeV. The only salient exceptions to this rule are the cases of 
A=82, 130, and 136. 

V. CONCLUSIONS 

In this work we have studied the 2vj3[3 decay matrix elements within a selfconsistent Skyrme Hartree-Fock calculation 
with pairing and deformation, and with residual ph and pp interactions treated in pnQRPA. The analysis has been 
focused on the study of the validity of SSD and LLSD. While the former has been tested on the 2v(3l3 decay of -'^""Mo, 
^^®Cd, and ^^^Te, the latter has been checked on all the measured 2vj3(3 emitters. Confirmation of that hypothesis 
has important consequences both experimentally and theoretically, since it could help to drastically simplify the 
experimental effort to get good estimates of 2v(3l3 half-lives on one hand, and to simplify the theoretical description 
of the intermediate nucleus on the other. 

With this aim, we have studied first the B{GT^) strength distributions of the intermediate 2^/3/? nuclei ^"^Tc, ^^^In, 
and ^^^I in the low range of excitation energies, comparing the theoretical results with the available experimental 
information from the direct decay and from charge- exchange reactions. The calculations reproduce fairly well the 
experimental distributions at low energy as well as the position and total strength of the GT resonances. 

The result of the analysis performed concludes that SSD is experimentally realized in the case of ^"''Mo because 
its experimental 2vj3[3 half-life is roughly accounted for by the measured B[GT) strengths connecting the ground 
states of the nuclei involved. On the contrary, the SSD half-lives are underestimated in the case of ^^^Cd and are 
overestimated in the case of ^^^Te. However, these results are still not conclusive since the uncertainties arising from 
the insufficient precision of the experimental measurements of ft-values, of 2v(3[3 half-lives, and especially because of 
the uncertainties of the B{GT) strength extracted from charge-exchange reactions. 

From the theoretical side, the conclusion of our results is that clear evidences for SSD and LLSD are not found 
within the present approach. This is at variance with the conclusion in Ref. p^ . where strong support of SSD was 
claimed. However, we would like to point out that our theoretical LLSD 2i>Pf3 matrix elements agree up to factors of 
the order of 2 with our total matrix elements. This finding is in qualitative agreement with the results obtained in 
previous studies [H, [l^. Furthermore, in Ref. [l5l] it is stated that SSD results agree quite well with the total matrix 
elements even though they find in some cases discrepancies up to a factor of 2. 

With only two exceptions (A=130, 136), our calculations reproduce the experimental 2i//?/3 half-lives when the whole 
energy range of excitation energies in the intermediate nucleus is considered. Within the considered nuclear model we 
find that important contributions to the 2ul3[3 decay matrix elements arise from relatively high excitation energies in 
several cases, not supporting dominance of single or even low-lying states. Further progress in both experimental and 
theoretical sides is still needed to clarify the importance of the contributions of the intermediate nucleus states to the 
2v(3(3 decay matrix elements. In the experimental case, this progress can be achieved by completing the experiments 
with high resolution charge-exchange reactions to more nuclei of interest and extending the measurements of GT 
strength distributions up to the resonance regions. For further progress in the field, measurements of the occupation 
numbers in particle-transfer reactions [STj are also of great importance. In the theoretical case improvements are 
needed in the treatment of the residual forces by using realistic interactions and improved parametrizations, especially 
in the particle-particle sector. The role of the axial-vector coupling strength should be also studied further because 
systematic indications in favor of strong quenching have been recently reported [38| . 

The single electron energy distributions of the outgoing electrons in the 2v(3l3 decay processes have been investigated 
under various assumptions (SSD, HSD), as well as from theoretical calculations using an exact treatment of the energy 
denominators of the perturbation theory, i.e. , without factorization of the nuclear structure and phase space integration 
calculations. A more comprehensive analysis of NEMO-3 data including the theoretical predictions presented in this 
paper could help to confirm or rule out the possible realizations of SSD and LLSD. 
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TABLE I: Experimental Gamow- Teller matrix elements M{GT ) and M{GT^) for the transitions connecting the T*" ground 
state in ^""Tc with the 0+ ground states in ^""Mo and ^°°Ru. We also show the matrix elements M'^^'^ [MeV"^] and the 
corresponding half-lives Ty2^ [y] obtained with aproximated (SSDl) and exact (SSD2) energy denominators. The experimental 
half-life is tII^I^ = 7.1 x 10^* y [H, and its corresponding matrix element is (MQy'''')exp = 0.241 MeV"^ While the half-lives 
are independent on gA, the matrix elements have been obtained using (;a=1.25. 

M{GT-) M{GT+) Mgy^'^ rf^f (SSDl) rf^f ^ (SSD2) 

0.57" 0.55* 0.19 1.1 xlO^^ 8.8 xlO^** 

0.65° 0.55'' 0.21 9.3 xlO^^ 7.2 xlO^'^ 

"from (3He,t) 
''from log /t^_ [32I 
'^from log /tEC U3 



TABLE II: Same as Table HI but for A = 116. In this case the 2i//3/3 half-life is T^J^I^ = 3 x 10^^ y [H, and its corresponding 
matrix element is {M'^^'^)^^-p = 0.127 MeV"^ 

M{GT-) M{GT+) M^r^'^ T^J^'^ (SSDl) T^f^'^ (SSD2) 

0.18" 0.51* 0.05 1.9 xlO^" 1.7 xlO^° 



0.51" 0.51' 0.14 2.4 xlO^^ 2.1 xlO^^ 
0.69^* 0.51' 0.19 1.3 xlO^^ 1.2 xlO^^ 



"from {^Re,t) 

''from log ft fl- [ 34l. used by [2I 
'^from (p,n)l23l_ 
''from log /tEC lI3| 



TABLE III: Same as TableH] but for A = 128. In this case the 2iy/3l3 half-hfe is T^J^'^ = 2.5 x 10^* y [H, and its corresponding 
matrix element is {Ml''^'^)^^p = 0.043 MeV"^ 

M(GT-) M{GT+) mI"^'^ T^f^ (SSDl) rf^;'^^ (SSD2) 

0.41" 0.10' 0.024 8.1 xlO^" 7.8 xlO^" 



0.33" 



0.10" 



0.019 



1.3 xlO-" 



1.2 xlO-" 



"from (p, n) [13 
''from log ftp- 
'^from log /tEC 
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TABLE IV: Theoretical (quenched) and experimental accumulated single Gamow- Teller strengths B{GT~) and B{GT^) up 
to the different measured excitation energies of the intermediate nucleus. The corresponding values of the double GT matrix 
element M^r^^ [MeV~^] up to the same excitation energies for ^°°Mo -^^""Ru are also given. 



E B(GT~) Yl B(GT+) Mpy^'^ 

E [MeV] th. exp." th. exp. th. exp 

1.4 0.42 0.46 0.73 - 0.12 

2.6 0.53 0.69 1.26 - 0.15 - 

'from (^'He.t) [3 



TABLE V: Same as in Table but for A = 116. 



J2 B(GT-') J2 B(GT^ 



M. 



2I//3/3 



E [MeV] th. exp." th. exp.' 



th. 



exp. 



0.7 
1.0 
2.2 
3.0 



"from (3He,t) | 
''from ((i,2He) 



0.29 0.03 

0.56 0.15 

0.56 0.32 

0.56 0.32 



0.50 
0.56 
0.69 
0.75 



0.33 
0.44 
<0.51 
<0.72 



0.03 
0.06 
0.06 
0.06 



0.05 
0.09 
<0.12 
<0.12 



TABLE VI; Same as in Table |IVl but for A = 128. 

J2 B(GT~) J2 B(GT+) Mgy^^ 
E [MeV] th. exp." th. exp. th. exp. 



0.58 


0.004 


0.22 


0.014 


0.000 


1.09 


0.094 


0.52 


0.021 


0.004 


1.61 


0.382 


0.64 


0.027 


0.012 


2.56 


0.468 


1.07 


0.187 


0.021 



"from (p, n) [3 
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FIG. 1: Accumulated B{GT~) strength for the transition ^""Mo^^™Tc (upper panel), accumulated B{GT+) strength for the 



transition ^""Ru 



"'°Tc (middle panel) and M^'^'''' [MeV"^] matrix element for the transition ^""Mo-^^'^Ru (lower panel) as a 
function of the excitation energy of the intermediate nucleus from a HF+BCS+QRPA calculations, together with experimental 
data from various sources. 
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FIG. 2: Same as in Fig. [T] but for the transition "®Cd-*"''Sn. 
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FIG. 3: Same as in Fig. [H but for the transition ^^^Te^^^^Xe 
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FIG. 6: Double-beta decay matrix elements as a function of the intermediate nucleus excitation energy considered in the calcu- 
lation. An experimental range is given by a shadow area where the upper and lower limits are deduced from the experimental 
half-lives 33], using =1.00 and qa =1.25, respectively. Dashed lines in = 48 and A — 76 correspond to calculations with 
an alternative choice of deformation (see text). 



